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Abstract 

Background: Tlie Hunt bumble bee {Bombus huntii Greene, Hymenoptera: Apidae) is a holometabolous, social 
insect important as a pollinator in natural and agricultural ecosystems in western North America. Bumble bees 
spend a significant amount of time foraging on a wide variety of flowering plants, and this activity exposes them to 
both plant toxins and pesticides, posing a threat to individual and colony survival. Little is known about what 
detoxification pathways are active in bumble bees, how the expression of detoxification genes changes across life 
stages, or how the number of detoxification genes expressed in B. huntii compares to other insects. 

Results: We found B. huntii expressed at least 584 genes associated with detoxification and stress responses. The 
expression levels of some of these genes, such as those encoding the cytochrome P450s, glutathione S-transferases 
(GSTs) and glycosidases, vary among different life stages to a greater extent than do other genes. We also found 
that the number of P450s, GSTs and esterase genes expressed by B. huntii is similar to the number of these genes 
found in the genomes of other bees, namely Bombus terrestris, Bombus impatiens, Apis mellifera and Megachile 
rotundata, but many fewer than are found in the fly Drosophila melanogaster. 

Conclusions: Bombus huntii has transcripts for a large number of detoxification and stress related proteins, 
including oxidation and reduction enzymes, conjugation enzymes, hydrolytic enzymes, ABG transporters, cadherins, 
and heat shock proteins. The diversity of genes expressed within some detoxification pathways varies among the 
life stages and castes, and we typically identified more genes in the adult females than in larvae, pupae, or adult 
males, for most pathways. Meanwhile, we found the numbers of detoxification and stress genes expressed by 
B. huntii to be more similar to other bees than to the fruit fly. The low number of detoxification genes, first noted 
in the honey bee, appears to be a common phenomenon among bees, and perhaps results from their symbiotic 
relationship with plants. Many flowering plants benefit from pollinators, and thus offer these insects rewards 
(such as nectar) rather than defensive plant toxins. 
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Background 

In addition to pesticide resistance, detoxification and 
stress responses are important adaptations that allow in- 
sects to overcome the chemical defenses of the plants 
and animals they feed on. Genes associated with these 
responses have been identified in many insects, includ- 
ing the mosquitoes Anopheles gamhiae [1] and Aedes 
aegypti [2], the fruit fly Drosophila melanogaster [3], the 
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honey bee Apis mellifera [4], and the red flour beetle 
Tribolium castaneum [5]. During pollination activities, 
bees are exposed to toxic substances in the environment, 
such as pesticides, phytochemicals, microbial toxins, 
pollutants and other xenobiotics [6-8], but genomic ana- 
lyses of the honey bee. A, mellifera, found fewer detoxifi- 
cation genes than are present in D, melanogaster and A, 
gamhiae [4,9]. To determine if this phenomenon is 
common among bees or unique to the honey bee, we 
evaluated the expression of detoxification and stress re- 
lated genes in a common western North American bum- 
ble bee, Bombus huntii (the Hunt bumble bee). Like A. 
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mellifera, B. huntii is a social insect with a diet based on 
pollen and nectar from a broad array of plants, so envir- 
onmental exposures to xenobiotics should be similar. In 
addition, bumble bee susceptibility to pesticides has 
been found to be similar to that of A. mellifera [10]. 

Generally, stress responses and the detoxification of 
xenobiotics includes three major and interrelated path- 
ways (Figure 1): oxidation-reduction, conjugation, and 
hydrolysis [11,12]. Oxidation-reduction enzymes include 
alcohol dehydrogenases, aldehyde dehydrogenases, cyto- 
chrome P450s, hydroxylases and peroxidases. Many 
cytochrome P450s are important enzymes for catalyzing 
oxidation-reduction reactions, and they are probably im- 
portant for pesticide detoxification in A, mellifera [13]; 
however, reactive oxygen species (ROS), such as hydro- 
gen peroxide, hydroxyl radicals, and superoxides, are 
typical by-products of these reactions. These ROS can 
be toxic in themselves, but are degraded by antioxidants, 
facilitated by other oxidation/reduction enzymes such as 
CuZn superoxide dismutases, catalases and peroxidases 
(Figure 1). In addition to the detoxification of insecti- 
cides and other xenobiotics, many of the oxidation/ 
reduction enzymes are also involved in the normal 
physiological functioning of insects. For example, cyto- 
chrome P450s are also involved in the production of 
pheromones, ecdysteroids and juvenile hormones 
[14-16], and many of the enzymes that breakdown ROS 
are associated with breaking down by-products from 
metabolism. 

The products of oxidation/reduction may be further 
degraded through conjugation. Conjugation is the 
process by which sugars, amino acids, or glutathione are 



transferred to a carrier compound for movement out of 
the cell. A host of conjugation enzymes catalyze these 
reactions [17] (Figure 1). Particularly important in in- 
secticide detoxification are members of the Glutathione 
S -transferases enzyme (GST) superfamily which catalyze 
the conjugation of oxidized lipids and exogenous toxins, 
such as pyrethroid and organophosphate pesticides 
[18-20] (Figure 1). GSTs are also involved in other func- 
tions, such as intracellular transport and the biosynthesis 
of hormones [21]. 

Detoxification is also carried out via hydrolysis, which 
transfers a hydrogen atom from water to the substrate 
compound, and these reactions are facilitated by hydro- 
lytic enzymes [12] (Figure 1). For example, carboxyl es- 
terases catalyze the hydrolysis of esters into an acid and 
an alcohol and are important in the degradation of pyre- 
throids and organophosphate insecticides [22]. In 
addition to these three pathways, other pathways may be 
invoked, such as those that involve ATP-binding cassette 
transporters (ABC transporters) [11], cadherins, or heat- 
shock proteins. 

Like all bumble bees, the life stages and castes of B, 
huntii vary considerably in morphology, behavior and 
activity, and so could be expected to differ in their ex- 
pression of detoxification and stress related genes. As 
with all Hymenoptera, bumble bees have a haplodiploid 
sex determination, where males are produced from 
unfertilized, haploid eggs and females from fertilized, 
diploid eggs [23]. Bumble bees live in annual, eusocial 
colonies founded by a single queen whose daughters 
form a female worker caste that provides foraging, brood 
care and nest maintenance [24]. A colony begins when a 
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Figure 1 Typical detoxification patiiways Icnown for insects. Solid boxes indicate chemical compounds from the environment or as a result 
of metabolic processes. Dashed-line boxes indicate the three detoxification pathways that this paper focuses on. 
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solitary overwintering queen emerges from hibernation, 
finds a suitable nesting site, provisions the nest with 
pollen and nectar, and commences laying eggs. The eggs 
hatch and larvae are fed nectar and pollen by the queen. 
In the fourth instar, the larvae begin to spin silk cocoons 
in which they pupate. The new adult bees emerge from 
their cells by chewing out of the cocoon, and these fe- 
male workers then take over the responsibilities of brood 
care and foraging as the nest grows. Near the end of the 
summer, the next generation of reproductive adults is 
produced and the colony senesces. The new queens fly 
from the colony to mate with males and then diapause 
as solitary individuals in a suitable location for the win- 
ter, whereas workers, males and the colony's original 
queen do not survive the winter. 

Our goal was to construct a transcriptome of the 
detoxification and stress response genes expressed in 
B. huntii. However, due to the complex life cycle of 
bumble bees, it is possible that not all the relevant genes 
are expressed in a single caste or life stage. To identify a 
fuller spectrum of genes expressed in B. huntii, messen- 
ger RNA (mRNA) was extracted from eggs, early instar 
larvae, late instar larvae, pupae, adult workers, adult 
males, an egglaying queen, and a diapausing queen. The 
corresponding cDNA was sequenced using pyrosequenc- 
ing, and genes associated with detoxification and stress 
response were identified. We also compared the number 
of cytochrome P450, GST and carboxylesterase genes 
found in B. huntii (acrosss all life stages) to those found 
in the fly D, melanogaster, and those found in other bees 
where annotated genomes are available, namely B. terres- 
triSy B, impatiens, A, mellifera, and M, rotundata, to 
evaluate whether the low number of detoxification genes 
found in A. mellifera [4] is unique to that bee, or more 
common among bees in general. 

Results 

An overview of detoxification and stress response genes 
in B. huntii 

We developed a cDNA database of 102,778 contigs and 
singletons using pyrosequencing of mRNA extracted 
from eggs, larvae, pupae, adult workers, adult males, an 
egglaying queen, and a diapausing queen of B, huntii. 
Gene expression data were examined two ways: the 
number of different genes detected within a class of 
detoxification genes (the gene diversity), and the number 
of transcripts of a particular gene or class of genes 
(the expression level). The total gene diversity within 
all the detoxification and stress response-related genes 
expressed in B. huntii was approximately 584 genes 
(Additional file 1: Table SI), after we removed sequences 
with high similarity to microorganisms and plants. 
Of the original detoxification and stress-related gene 
sequences we identified in the samples, 13.8% were 



bacterial, 12.4% were similar to known plant sequences 
and 3.6% were similar to known fungal sequences. These 
were all removed. Among the remaining genes, we identi- 
fied genes associated with oxidation-reduction, conjuga- 
tion, and hydrolytic enzymes, as well as some other 
detoxification and stress responses (Additional file 1: 
Table SI, Tables 1 and 2). The genes for enzymes in each of 
these four groups were expressed differently in the various 
life stages and castes, with some expressed in multiple stages 
and others expressed in only one or two (Tables 1 and 2). 

We found 202 genes associated with oxioreductases. 
Most of these were expressed across the majority of the 
life stages, but adult females had the greatest diversity of 
expressed gene types in this category (Table 2). 

Of the 109 conjugation-enzyme related genes we iden- 
tified, GSTs, glycosyltransferases, methyltransferases and 
aminotransferases (transaminases) were expressed in all 
life stages. The genes coding for acetyltransferases, acyl- 
transferases, CoA transferases, formyltransferases, and 
phosphotransferases were also expressed in different life 
stages, but abundantly in adult females (Figures 2 
and 3). Genes encoding sulfo transferases were expressed 
only in the diapausing queen (Figures 2 and 3). 

We detected a total of 143 genes for hydrolytic en- 
zymes, across all castes and life stages (Tables 1 and 2), 
but the greatest diversity of expressed gene type was 
again found in the adult females (Table 2). 

In addition to the genes in these three major groups, 
we found 130 other genes encoding proteins that poten- 
tially have a detoxification or stress response function, 
such as ABC transporters, cadherins, heat shock pro- 
teins, isomerases and lyases (Additional file 1: Table SI). 
Late instar larvae, pupae and adult males had the lowest 
diversity of these genes (Tables 1 and 2). 

Gene expression among the different life stages of B, 
huntii 

The goal of this project was to survey detoxification and 
stress related genes in B. huntii. To do this, we used about 
1000 eggs, 20-50 bees for the immature stages (see 
Methods), and one or two bees for the adult stages, and 
then pooled samples within each life stage for the 454 se- 
quencing. Our sampling, which was not truly replicated 
within each life stage, provides some preliminary informa- 
tion on the expression levels among the different life 
stages and castes, but no statistical comparisons could be 
made. Among the immatures, the oxidoreductases were 
well expressed in the larvae, especially the late instar lar- 
vae (Figure 2). Within the conjugation enzyme group, the 
GSTs were highly expressed in the immatures, especially 
the pupae, but not the eggs (Figure 2). Genes associated 
with hydrolytic enzymes were also highly expressed in the 
immatures, especially the nitrilases, and again, particularly 
in the late instar larvae sample (Figure 2). 
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Table 1 Putative functions for the detoxification genes expressed in different life stages of Bombus huntii, based on 
the function assigned when the gene was first isolated 



Enzymes 



Putative functions 



Taxon of identified gene Reference Bombus 

huntii life 
stage/caste^ 



Oxidation and reduction enzymes 



Aldehyde 
dehydrogenases 

Alcohol 

dehydrogenases 
Catalases 



Acetaldehyde detoxification 



Interconversion of alcohols and aldehydes or ketones with the 
reduction of NADH 

Catalyze H2O2 to H2O and O2 



Drosophilo melonogoster [25] 
D. melonogoster [25] 
Anopheles gombloe [26] 



II, III, V, VI, VII 



I, II, III, VI, VII, 

viii 



Cytochrome P450s Detoxify xenobiotics by transferring O2 to a substrate and 
producing H2O2, OH", O2 



Dehalogenases 
Hydroxylases 

Oxidoreductases 

Peroxidases 

Superoxide 
dismutases (SOD) 

Thioredoxins/ 
glutaredoxins 

Conjugation enzymes 



Remove halogen from haloacid compounds 

Degrades toxic organic compounds and pesticides by 
introducing -OH 

Detoxify gramine 

Catalyze H2O2 and organic hydroperoxides (OHPs) 
Catalyze the dismutation of O2 into O2 and H2O2 

Detoxification of ROS 



D. melonogoster, A. gombloe, [1,16,27,28] all 

Tyrio jocoboeoe, Culex 

quinquefosclotus 

SInorhlzoblum mellloti [29] 

Chrysopogon zlzonloldes [30] 

SItoblon ovenoe [31] 

A. gombloe [32] 

A. melllfero [33,34] 

Socchoromyces cerevlsloe [35] 



Acetyltransferases 

Acyltransferase 

CoA transferases 
Formyltransferases 

Glutathione S- 
transferases (GST) 

Glycosyltransferases 

Methyltransferases 

Phosphotransferases 

Sulfotransferases 

Transaminases or 
aminotransferases 

Hydrolytic enzymes 

Acetylcholine and 
carboxyl esterases 

Acid/Alkaline 
phosphatases 

Amidases 
Aminopeptidases 
Cyclohydrolase 



Detoxify sulfadiazine and aromatic chemicals by transferring 
acetyl groups from acetyl-CoA to arylamines 

Catalyze carboxylic acid group, such as benzoic, isovaleric, or 
acetylsalicylic acids 

Acetate detoxification 

Formaldehyde detoxification 

Detoxify endogenous compounds (peroxidised lipids) and 
exogenous toxins (pyrethroid and organophosphate) by 
catalyzing reduced glutathione 

Glycosyltransfer to 0H-, NH2-, SH-, or COOH for detoxifying 
insecticides 

Catalyze methylation reactions using S-adenosyl-L-methionine as 
a substrate for detoxification 

Detoxify insecticides such as DDT 

Detoxify insecticides such as DDT 

Detoxify 3-hydroxykynurenine and glyoxylate by catalyzing 
amino acid and a-keto acid 



Detoxify pyrethroids and organophosphates by hydrolysis 
of acyl group or ester bonds 

Detoxify endotoxin and plant toxin 



Detoxify toxic amides or esters, or organophosphorus or 
carbamate pesticides 

Detoxify Bt Cry toxin, mycotoxin, organophosphonates, 
pyrethroid esters, microbial or botanical pesticides 

Hydrolysis of L-phenylalanine 



Leglonello pneumophllo. Homo [36,37] 
soplens 

H. soplens [38] 

Aspergillus nidulons [39] 

Burkholderlo fungorum [40] 

Tenebrlo molltor, Bombus [18,41] 
Ignites, Locusto migrotorio 
monllensis 

Bombyx morl [17,42] 

H. soplens [43] 

Triotomo Infestons [44] 

T. Infestons [44] 
Aedes oegypti, D. melonogoster [45] 



III, IV, VI, VII 

ii-vii 

all 
all 
all 

i, ii, vi, vii, viii 



I, II, III, V, VI, VII, 

viii 



I, II, IV, VI, VII, 

viii 



I, II, III, VI, VII, 

viii 



all 



all 

i, ii, vi, vii, viii 
iii, vii 
ii-viii 



Lucilio cuprino, D. [22] all 
melonogoster 

Donio rerlo, Agrlollmox ogrestis [46,47] ii, iii, vi, vii 

Acromyrmex octosplnosus [48] 

B. morl, Pseudopluslo [42,48,49] all 
includens, E. coll 

Drosophilo spp. [50] vi 



II, III, IV, VI, VII, 

viii 
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Table 1 Putative functions for the detoxification genes expressed in different life stages of Bombus huntii, based on 
the function assigned when the gene was first isolated (Continued) 



Glycosidases 


Hydrolyze sugar-containing compounds such as tannic acid and 
gallic acid 


Manduca sexto 


[48] 


i, ii, iii, v, vi, vii, 
viii 


Glycosylases 


Detoxify cytotoxic and cytostatic substances such as 5-methyl 
group of a thymine residue, biodegrade saponins 


Phaeosphaeria nodorum, P. 
ovenorio f. sp. triticeo 


[51] 


i, ii, vi, vii, viii 


Nitrilases 


Detoxify HCN 


Caenorhabditis elegans 


[52] 


i-vii 


Phosphodiesterase 


Break phosphodiester bond for hydrolysis of organophosphate 
insecticides 


E. coll 


[53] 


ii, iii, iv, vi, vii 


Phosphohydrolases 


Detoxify cypermethrin and bifenthrin 


Tribolium castaneum 


[54] 


ii, iii-viii 


Others with possible stress related functions 








ABC transporters 


Export conjugated toxins out of the cell 


Drosophila spp. 


[11] 


vii, viii 


Cadherins 


Resistance to Bt toxin (Cry 1 Ac) 


Helicoverpa armigera 


[55] 


i, ii, vi, vii, viii 


Heat shock proteins 


Involved in folding and unfolding of proteins, detoxification of 
pesticides and heavy metals 


Chironomus tentons 


[56] 


all 


Isomerases 


Detoxify organophosphates by structural rearrangement of 
isomers 


E. coli 


[16] 


all 


Lyases 


Break various chemical bonds other than hydrolysis and oxidation 
for detoxification 


Zygoeno tronsolpino 


[57] 


all 



i-egg, ii-early instar larvae, iii-late instar larvae, iv-pupae, v-adult male, vi-adult worker, vii-diapausing queen, viii-egglaying queen, all-all stages. 



Genes associated with oxidation and reduction were 
well expressed in all the adult stages. The egglaying 
queen had quite high expression levels of these genes, 
although not for all the genes (Figure 3). As with the im- 
matures, GSTs were the most well expressed genes in 
the conjugation enzyme group. In contrast to the high 
GST expression in the immatures, the glycosidases in 
adults were the most expressed of the hydrolytic en- 
zymes, particularly in the adult workers and the egglay- 
ing queen (Figure 3). Heat shock protein genes were 
well expressed in the eggs, the diapausing queen and the 
egglaying queen, while expression levels in the males 
was fairly low (Figure 4). 

Comparison of detoxification genes in B. huntii with other 
insects 

Gene diversity for detoxification and stress related genes 
in B. huntii was similar to, or lower than, that for the 



other bees {B. terrestris, B. impatiens, A. mellifera and 
M, rotundata), and generally a much lower gene diver- 
sity than D, melanogaster (Table 3). Drosophila melano- 
gaster has many more P450s, acyltransferases, GSTs, and 
most of the hydrolases, than any of the bees. One excep- 
tion is the nitrilases, for which we found a high diversity 
in B, huntii. The high nitrilase gene diversity does not 
occur in the other bees (Table 3). Our data was based on 
a transcriptome and may underrepresent the number of 
genes, relative to the number identified in the genomes 
of the other species included in this analysis. 

We also compared some subgroups within the cyto- 
chrome P450s and GSTs. Within the cytochrome P450s, 
the CYP2 and mitochondrial P450 enzymes tend to have 
endogenous functions, while the CYP4 and CYP6 en- 
zymes have detoxification functions. We found a similar 
number of CYP6 related genes among the six insects, 
however a greatly reduced number of CYP4 related 



Table 2 The number of detoxification genes identified, as expressed in different life stages of Bombus huntii 


Stages 


Oxidation-reduction genes 


Conjugation genes 


Hydrolysis genes 


Other detoxification genes 


Total 


Egg 


36 


14 


22 


21 


93 


Early instars 


49 


20 


33 


23 


125 


Late instars 


64 


18 


35 


12 


129 


Pupa 


31 


9 


13 


7 


60 


Adult male 


36 


11 


23 


13 


83 


Adult worker 


87 


40 


59 


60 


246 


Diapausing queen 


93 


40 


54 


60 


247 


Egglaying queen 


68 


30 


36 


47 


181 


All stages 


202 


109 


143 


130 


584 
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Figure 2 The normalized expression levels of possible detoxification genes in B. huntii immatures, for genes in the (A) oxidation: 
reduction, (B) conjugation, and (C) hydrolysis groups. All data were normalized to a million mapped reads per sample (1 RPKM = number of 
reads per million mapped reads of genes). 



genes in the bees, relative to D. melanogaster (Table 3). 
The GSTs are composed of several sub-groups, but the 
sigma-GSTs and epsilon-GSTs comprise the main sub- 
groups with detoxification functions. We found very few 
sigma-GSTs occurred in any of these insects; however, 
fewer epsilon-GSTs occurred in the bees than D, mela- 
nogaster (Table 3). 



Discussion 

Female adults (workers and queens) had the greatest di- 
versity and gene expression levels among the detoxifica- 
tion and stress related genes we identified in B, huntii. 
Female adults expressed high levels of putative detoxifi- 
cation genes, while adult males and pupae had relatively 
lower expression levels of these genes. The low diversity 
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Figure 3 The normalized expression levels of possible detoxification genes in B. huntii adults, for genes in the (A) oxidation: reduction, 
(B) conjugation, and (C) hydrolysis groups. All data were normalized to a million mapped reads per sample (1 RPKM = number of reads per 
million mapped reads of genes). 



and expression levels in males may in part be associated 
with haploidy, as was demonstrated in the stingless bee 
Melipona quadrifasciata [58]. Behaviors associated with 
the different bumble bee castes may also affect gene ex- 
pression. Male activities are all related to feeding them- 
selves and mating; whereas, both workers and queens 
have more complex behaviors, often with high energy 
demands (e.g. foraging, brood care, and egglaying). 



Increases in foraging and feeding activity directly in- 
crease the risk of exposure to environmental pollutants; 
furthermore, increased energy demands have been 
shown to increase consumption rates and the formation 
of metabolic by-products in A, mellifera [59]. An actively 
laying queen may not forage after the first batch of 
brood matures, but she does have high metabolic de- 
mands for egglaying, policing the nest to keep workers 



Xu et al. BMC Genomics 2013, 14:874 
httpy/www.biomedcentral.com/l 471 -21 64/1 4/874 



Page 8 of 1 3 




Figure 4 Normalized expression levels of heat shock protein 
genes in B. huntii across all the life stages and castes. All 

data were normalized to a million mapped reads per sample 
(1 RPKM = number of reads per million mapped reads of genes). 



from becoming egg layers, and incubating the brood. 
And, indeed, the diversity and expression of potential 
detoxification genes is relatively high in queens and 
workers. 

Additionally, detoxification activity could be affected 
by variations in hormone levels and the complex mor- 
phological changes that occur in holometabolous insects. 
Changes of hormone levels during insect development 
are well quantified in Drosophila, with the steroid hor- 
mone ecdysone showing peaks during the transition 
from larva to pupa [60]. Metabolic activity also varies 
between life stages; most notably it declines during dia- 
pause, a dormant state during which respiration is very 
low. Thus, the fact that detoxification genes are more 
highly expressed in an egglaying than a diapausing 
queen is not entirely surprising. Non-diapausing adults 
are more active than pupae, but pupae undergo a major 
metamorphosis, so the relative expression of detoxifica- 
tion and stress related genes is more difficult to predict, 
but is more likely to be the result of metabolic processes 
and not a result of xenobiotic exposure. 

Some detoxification mechanisms may be required 
more during certain developmental stages, and some 
genes that we classified as potential detoxification genes 
may serve non-detoxification functions. For example, 
glycosidases are more highly expressed in queens and 
workers than in other stages. Glucosidase catalyzes the 
cleavage of individual glucosyl residues from various 
glyco-conjugates, a process involved in the breakdown 
of sugar-containing compounds during the digestion of 
pollen and nectar. GSTs also break down a variety of 
compounds produced during metabolism and digestion, 
such as partially digested lipids, in addition to their in- 
volvement in the detoxification of environmental toxins, 
plant allelochemicals, and organochlorine and organo- 



phosphate insecticides [18-20]. GST genes are expressed 
at high levels in B. huntii adults, pupae, and larvae, but 
at low levels in eggs, which may be related to either dif- 
ferences in metabolism or exposure to environmental 
toxins. For example, late instar bee larvae accumulate 
high levels of lipids [61,62] and the high expression of 
GSTs in pupae may be related to the metabolism of 
these lipids during pupation. Similar results have been 
reported in other insects such as A, aegypti, Lucilia 
cuprina, and Tenebrio molitor, where the enzyme activ- 
ities of GSTs were particularly high in the pupal stage 
[63]. In addition, some cytochrome P450s are involved 
in lipid metabolism, perhaps explaining why late instar 
larvae had higher expression levels of cytochrome P450s 
than early instar larvae. 

Our gene analysis was based on a transcriptome, and 
as such, might underrepresent the number of detoxifica- 
tion genes found in the genome. However, we found 
similar numbers of genes to those found in the genomes 
of other bees. In addition, we used a large number of in- 
dividuals for our transcriptome, and from across all the 
life stages of the bee. Thus, our transcriptome is prob- 
ably a good library of the detoxification and stress re- 
lated genes found in B, huntii. The fact that our bees 
were produced in culture on sugar syrup, rather than 
flower nectar, may have reduced the expression of some 
detoxification genes, such as has been found with honey 
bees [64]. However, we also fed our bees pollen collected 
by honey bees, and this likely provided a source of plant 
based phenolics and perhaps even traces of pesticides. 

Due to the importance of detoxification capabilities 
during a bumble bee s life, one might expect B, huntii to 
have a genome rich in detoxification genes; however, like 
other bees such as A. mellifera, B, terrestris and B. impa- 
tiens, we found B, huntii to have relatively few P450s, 
GSTs and carboxylesterases/esterases, as compared to D, 
melanogaster. For example, the number of P450 CYP4 
genes is quite high in D, melanogaster, relative to what 
we found in B, huntii, Bombus huntii also has fewer 
epsilon-GSTs genes in comparison to D, melanogaster, 
although we found a slightly higher number of sigma- 
GSTs in B, huntii and A, mellifera, Epsilon-GSTs are 
known to be involved in DDT (dichlorodiphenyltrichloro- 
ethane) resistance and the detoxification of xenobiotics. 
The high number of epsilon-GSTs in D, melanogaster 
must reflect either a greater evolutionary diversification of 
these genes in the fruit fly, a loss of some genes in the 
bees, or a combination of both. In any case, this could 
arise if flies tend to have a greater exposure to environ- 
mental toxins than do the bees. Sigma-GSTs are associ- 
ated with the detoxification of Upid peroxidation 
products, and provide protection from oxidants pro- 
duced by the aerobic metabolism of honey, pollen and 
nectar in the bees [9]. Thus, the slightly greater number 
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Table 3 Comparison of the number of detoxification genes identified in Bombus huntii, in other bees, and in 
Drosophila melanogaster 



Putative function 



Bombus 
huntii 



Bombus 
terrestris 



Bombus 
impatiens 



Apis 

mellifera 



Megacliile 
rotundata 



Drosophila 
melanogaster 



Oxidation and reduction enzymes 



Aldehyde dehydrogenases 

Catalases 

Cytochrome P450 

(CYP4) 

(CYP6) 

Glutaredoxins 
Peroxidases 

Superoxide dismutases 
Conjugation enzymes 

Acyltransferases 

Glutathione S-transferases 

(Sigma-GSTs) 

(Epsilon-GSTs) 

Sulfotransferases 

UDP- 

glucuronosyltransferases 
Hydrolytic enzymes 

Alkaline phosphatases 
Amidases 
Aminopeptidases 
Amylases 

Carboxylesterases/Esterases 

Glucosidases 

Glucuronidases 

Glycosylases 

Nitrilases 



5 
2 

44 

(2) 

(18) 

2 

16 

2 

13 

11 

(3) 

(0) 

2 

2 



1 

2 

12 

2 

23 

8 

1 

3 

25 



5 
1 

50 

(6) 

(22) 

3 

9 

4 

28 

14 
(0/?) 
(0/?) 
17 



3 
4 
19 
1 

17 

7 

1 

2 

1 



5 

2 

49 

(5) 

(22) 

7 

9 



30 
15 

(0/?) 
(0/?) 
16 



26 
1 

22 

8 

1 

5 

2 



46 

(5) 

(28) 

4 

17 

5 

32 
18 
(4) 
(0) 
17 
2 



3 
4 
24 
2 

24 

18 

1 

7 

3 



18 
5 

52 
(6) 
(19) 
3 



9 

(1) 
(0) 
5 
2 



2 
7 

27 
2 

22 

5 

2 



23 
3 

85 
(22) 
(22) 
7 

20 
6 

55 

42 

(1) 

(10) 

15 

7 



16 

12 

41 

15 

54 

29 

4 

5 

2 



Note: The numbers of B. huntii and M. rotundata genes are transcriptome data from 454 pyrosequencing, except for P450 of M. rotundata which Is genomic data. 
The numbers of B. terrestris, B. impatiens, A. mellifera, and D. melanogaster genes are from genomic data from the NCBI GenBank database (http://www.ncbi.nlm.nih.gov/ 
genome, accessed on September 25, 2013), and the numbers in brackets below P450 or GSTs are subgroups of P450 or GSTs. 



of sigma-GSTs in bees also appears to reflect their eco- 
logical niche. 

We found here that the low number of detoxification 
genes found in the A. mellifera genome [4] is a 
phenomenon common to many bees. We hypothesize 
that this phenomenon may have arose as a result of the 
symbiotic relationship between bees and flowering 
plants. Flowering plants often produce rewards (such as 
nectar) to attract bees and other pollinators. Further- 
more, at least some plants have lower levels of plant de- 
fensive compounds in the pollen and nectar [8], and 
adding plant alkaloids to the nectar reduces pollinator 
activity on those flowers [65]. Thus, the detoxification 
abilities of bees may be less than the flies (and likely 
other insects as well) due to a lower level of exposure to 
plant defensive compounds, compounds that plants 



produce to defend themselves against herbivores not 
pollinators. 

Conclusions 

Using transcriptome analysis of all life stages, we found 
the Hunt bumble bee, B, huntii, to have the genetic po- 
tential to produce a large number of detoxification and 
stress related proteins, including oxidation and reduc- 
tion enzymes, conjugation enzymes, hydrolytic enzymes, 
ABC transporters, cadherins, and heat shock proteins. 
The number of genes in these pathways was fewer than 
found in flies, such as D. melanogaster, and slightly 
lower than that found in the bumble bees B. terrestris 
and B. impatiens, the honey bee A. mellifera, and the 
solitary bee M rotundata. However, a transcriptome 
may underestimate gene diversity, as compared to 
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studies based on a genome. We also found that, in gen- 
eral, low levels of detoxification and stress related genes 
are expressed in pupae, adult males and larvae than in 
adult females. Workers and queens express high levels 
of P450s and glycosidases. 

Methods 

Source of B. huntii 

Eight different stages of B. huntii were used in this ana- 
lysis: eggs, early (2nd and 3rd) instar larvae, late (4th) 
instar larvae, pupae, adult workers, adult males, a dia- 
pausing queen, and an egglaying queen. All stages were 
collected from a nest cultured in the lab at the USDA- 
ARS Pollinating Insect Research Unit in Logan, UT, ex- 
cept for the diapausing queen, which was a sister of the 
egglaying queen and had been held in cold storage at 
4°C for three months prior to collection for sequencing. 
The bees were reared according to Strange [66] and 
were started from queens that were raised and mated in 
the laboratory. The colony was fed on a diet of pollen 
collected from honey bee colonies and a 1:1:2 glu- 
cose:fructose:sucrose syrup solution. The eggs, larvae 
and pupae were removed from the colony and killed 
directly by immersion in ^NAlater solution (Life 
Technologies, NY, USA), whereas the adult bees were 
first killed by immersion in liquid nitrogen and were 
then placed in vials of RNAlater solution. All bee tis- 
sues were submerged in approximately 5 volumes of 
RNAlater solution and kept at 4°C overnight to per- 
meate the cells for stabilizing the RNA, the samples 
were then stored (about one month) at -80°C until 
processed. 

Preparation of RNA 

For RNA isolation, about 1000 eggs, 50 2nd-3rd instar 
larvae, 20 4th instar larvae, 20 pupae, two adult males, 
two adult workers, one egglaying queen and one 
diapausing queen were removed from RNAlater and 
washed twice with nuclease free water and then trans- 
ferred to a mortar and ground in liquid nitrogen to a 
fine powder. The total RNA was extracted by resuspend- 
ing the ground powder into 20 ml extraction buffer 
[100 mM NaCl, 2% SDS, 50 mM Tris-HCl (pH 9.0), 
10 mM EDTA (pH 8.0)] and 20 ml phenol-chloroform- 
isoamyl alcohol (lAA) (49.5:49.5:1, by volume) in a 
50 ml centrifuge tube. The solution was mixed and then 
centrifuged at 8,000 rpm for 20 min at 4°C. The aqueous 
phase was removed and placed in a clean centrifuge 
tube and an equal volume of phenol-chloroform-IAA 
(49.5:49.5:1) was added. The mixture was shaken and 
then centrifuged at 8,000 rpm for 20 min at 4°C. This 
organic extraction was repeated two more times. The 
RNA was precipitated with a 1/10 volume of 3 M so- 
dium acetate (pH 5.2) and 2.5 volumes of 95% ethanol. 



The RNA pellet was washed with 70% ethanol, dried for 
5 min, and resuspended in 400 [A RNase-free water con- 
taining 1% diethylpyrocarbonate (DEPC). We repeated 
total RNA extraction once for each treatment. 

Poly(A) + mRNAs were purified with an oligo(dT)-cel- 
lulose column through the binding, washing and elution 
steps. First, 1 ml of total RNA solution (1 mg of RNA 
for each sample, except -200 [ig for eggs, this was half 
of the RNA from each extraction) was heated at 65°C for 
5 min, then cooled on ice for 5 min, and 200 [A sample 
buffer [3 M NaCl, 10 mM Tris (pH 7.5), 1 mM EDTA] 
was added. For the binding step, 8.8 ml of binding solu- 
tion [0.5 M NaCl, 10 mM Tris (pH 7.5), 0.5% SDS, 
0.1 mM EDTA] was added to 1.2 ml RNA sample, agi- 
tated for 30 min and then briefly centrifuged to remove 
the supernatant; all steps were repeated twice more. For 
the washing step, 10 ml of high salt buffer [0.5 M NaCl, 
10 mM Tris (pH 7.5), 1 mM EDTA] were added to the 
oligo(dT) -cellulose, which was then mixed by rotating 
2 min, followed by a brief centrifugation to remove the 
supernatant. The oligo(dT) was then suspended in 
10 ml of high salt buffer and transferred to a 20 
ml-column (Bio-Rad, Hercules, CA, USA), washed with 
the high salt buffer twice, then washed another time 
with a low salt buffer [0.5 M NaCl, 10 mM Tris 
(pH 7.5), 1 mM EDTA]. Pre-warmed (65°C) elution buf- 
fer [10 mM Tris (pH 7.5), 1 mM EDTA] (3 ml) was 
added to the top of the oligo(dT)-cellulose for a third 
time, the suspension was collected, and mRNA was 
precipitated by adding 50 [A of glycogen solution 
(20 mg/ml), 1/10 volume of 3 M NaAc (pH 5.2), 
7.5 ml of 100% chilled ethanol, and stored overnight 
at -20°C before being centrifuged. The mRNA pellet 
was washed with 70% ethanol and dried for 10 min, 
and then dissolved in 80 [A of RNase-free water (1% 
DEPC). In addition, the mRNA samples from eggs, 
larvae, pupae and adult males were amplified using 
the MessageAmp III RNA amplification kit (Ambion, 
Austin, TX, USA), creating a sample that was cRNA. 

cDNA library preparation for 454 sequencing 

For each sample, a cDNA library was prepared with 
mRNA or cRNA using a cDNA rapid library prepar- 
ation kit (Roche, Branford, CT, USA) according to 
the manufacturers instructions, with minor changes. 
Briefly stated, 18 [A of mRNA or cRNA (around 
500 ng) were fragmented using fragmentation solution 
[0.1 M ZnCl2, 0.1 M Tris-HCl (pH 7.0)], followed by 
vortex mixing and a brief centrifugation, and then 
heated at 70°C for 30 s. The reaction was stopped by 
chilling on ice and adding 2 \A of 0.5 M EDTA 
(pH 8.0) and 28 \A of 10 mM Tris HCl (pH 7.5). The 
mRNA fragments were purified using 80 [A of 
RNA Clean reagent, containing SPRI beads (Beckman 
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Coulter, Beverly, MA, USA), and 19 \i\ of 10 mM 
Tris HCl (pH 7.5). The beads were removed with 
centrifugation, and the supernatant containing the 
RNA was added to a new 200 [A tube. 

The first-strand cDNA was synthesized by adding 8 [A 
5x RT-buffer AMV, 4 (il 0.1 M DTT, 4 [A 10 mM dNTP, 
1 [A protector RNase inhibitor (25 U/(il), 2 [A AMV RT 
(25 U/[A) to the clean, fragmented RNA, gently mixing, 
then incubating at 25°C for 10 min, followed by 42°C for 
60 min. The second strand cDNA was synthesized by 
mixing in 30 (il 5x second strand synthesis buffer, 1.5 [A 
10 mM dNTPs, 6.5 \A 2^"^ strand enzyme and 72 \A 
double distilled water before incubating at 16°C for 2 h, 
then adding 20 [A T4 DNA polymerase, incubating at 
16°C for 5 min, and finally adding 17 [A of 0.2 M EDTA 
(pH 8.0) to stop the reaction. The double stranded 
cDNA was purified using AMPure beads (Beckman 
Coulter, Beverly, MA, USA), and the cDNA was then 
dissolved in 16 [A of 10 mM Tris-HCl (pH 7.5). The 
cDNA was further purified using gel purification to iso- 
late fragments of 500-800 bp. 

To repair fragment ends, 9 [A of end repair mix (2.5 [A 
RL (Rapid library) lOx buffer, 2.5 [A RL ATP, 1.0 [A RL 
dNTP, 1.0 \A RL T4 polymerase, 1.0 \A RL PNK and 
1.0 [A RL Taq polymerase), from a cDNA RL preparation 
kit (Roche, Branford, CT, USA), were added to the 
cDNA, incubated at 25°C for 20 min, 72°C for 20 min, 
and then held at 4°C. The adaptor ligation was com- 
pleted by adding 1 [A of RL adaptor and 1 [A of RL ligase 
to the reaction tube and incubating at 25°C for 10 min. 
The small fragments (less than 100 bp) were removed 
using AMPure beads, and the supernatant contained the 
cDNA library. The cDNA libraries were then amplified 
by running emulsion-PCR and sequence analysis per- 
formed on a Roche GS-FLX system at the Center for In- 
tegrated BioSystems (CIB), Utah State University, Logan, 
Utah (http://biosystems.usu.edu/). 

Sequence assembly, annotation and detoxification gene 
identification 

The 454 sequence outputs were aligned and assembled 
de novo using CLC Genomics Workbench. The contigs 
and singletons obtained from de novo assemblies were 
BLAST-searched against the GenBank database at the 
National Center for Biotechnology Information (NCBI) 
(http://www.ncbi.nlm.nih.gov/BLAST/) in the iNquiry 
Bioinformatics Portal. Those similarities maintaining E- 
values of less than or equal to 0.001 (more than 93% of 
matches were based on E < e'^^) were treated as signifi- 
cant matches and were selected as the annotation of B, 
huntii unigenes. The detoxification genes were identified 
by comparison with detoxification genes found in A, 
melliferuy D, melanogaster and other organisms. 



The quality measurement of RNA, cDNA and sequences 
assembly 

The quality and integrity of total RNA, mRNA and 
cDNA is very important for obtaining high quality tran- 
scriptome sequences. The concentration of total RNA 
was measured using a NanoDrop 2000 Spectrophotometer 
(Thermo Scientific, Inc. Wilmington, DE, USA), and the 
quality and integrity of total RNA was examined by elec- 
trophoresis. We only used total RNA samples with an 
A260/A280 ratio of 2.0 to 2.2 and two typical rRNA bands. 

The mRNA concentrations were measured after the 
first purification from total RNA, and then again after 
fragmentation for cDNA preparation, using a spectro- 
photometer (Turner BioSystems, Inc. Sunnyvale, CA 
USA) with ribogreen RNA reagent (Invitrogen, Foster, 
CA, USA). The quality and integrity of the mRNA was 
examined using an Agilent 2100 Bioanalyzer with RNA 
6000 Pico kit (Agilent Technologies, Santa Clara, CA, 
USA). All samples used were of high quality and integrity, 
as determined by mRNA fluorescence figure with typical 
shape of broad peak and without two ribosomal RNA con- 
tamination peaks. The quality of the fragmented mRNA 
were determined by running 1 \A of the fragmented 
mRNA and 1 \A of non-fragmented mRNA on an RNA 
6000 Pico Chip on the Agilent 2100 Bioanalyzer. All frag- 
mented samples showed lengths of around 800 bp. 

The quality of the cDNA library was determined by 
the Center for Integrated BioSystems (CIB), Utah State 
University using a high sensitivity DNA assay on an 
Agilent Bioanalyzer. All samples displayed a broad shape 
of peak from 600 bp to 1200 bp with a relatively higher 
peak at approximately 800 bp. The total read count for 
the 454 sequence after assembly was 837,010. The aver- 
age read length was 425 bp with a total read length of 
355,789,178 bp. 

Gene diversity and expression levels for detoxification 
and stress related genes 

The diversities of detoxification genes were determined 
by identifying the number of genes in the particular en- 
zyme group using a BLAST search against the GenBank 
database at NCBI. Assembled contigs from B, huntii that 
differed from each other in sequence, but matched the 
same gene in GenBank were considered to be different 
regions of the same gene if the contigs were each shorter 
than half the sequence length of the GenBank gene, 
otherwise they were considered to be different genes. 

The expression levels of individual detoxification genes 
were estimated using RNA-seq [67,68] as follows: 



Expression level = (Total reads of the gene) 

/ (Total mapped reads in millions for the sample 
X gene length in Kb) 
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The cumulative expression level for a group of genes 
was calculated as follows: 

Expression level = (Total reads within the group of genes) 

/ (Total mapped reads in millions for the sample 
X average length of the group of genes in Kb) 
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Additional file 1: Table SI. Detoxification and stress response related 
genes identified in Bombus huntii. 
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